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Abstract 

In spite of precautions to avoid the harmful effects of extreme events, we experience 
recurrently phenomena that overcome the preventive barriers. These barriers usually 
increase drastically right after the occurrence of such extreme events, but steadily 
decay in their absence. In this paper we consider a simple model that mimics the 
evolution of the protection barriers to study the efficiency of the system's reaction 
to extreme events and how it changes our perception of the sequence of extreme 
events itself. We obtain that the usual method of fighting extreme events introduces 
a periodicity in their occurrence and is generally less efficient than the use of a 
constant barrier. On the other hand, it shows a good adaptation to the presence of 
slow non-stationarities. 
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1 Introduction 



One important motivation for the unified study of extreme events is the con- 
centration of the destructive power of different systems in some few rare events. 
Remarkable examples are earthquakes, extreme weather conditions, epileptic 
seizures, heart attacks, stock markets crashes, etc... (Q). Extreme events are 
usually generated by complex dynamics that involves the coupling between 
many different dimensions and scales (J2|). However, the characterization of 
extreme events is usually done in a single scientifically or socially relevant 
observable, like the magnitude of earthquakes, number of days of drought, or 
the highest wind speed in storms. From the point of view of dynamical sys- 
tems, which is assumed throughout this paper, we say that these observables 
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are obtained applying an observation function to the full phase space of the 
system. 

In this paper we assume a general perspective for the influence of feedback 
reactions to extreme events. These reactions may be consequence of human 
activities or of natural feedback loops present in the system (J3j). For concrete- 
ness we motivate the problem and we associate our mathematical model with 
the case of human reactions to floods in rivers (0), which is a representative 
example of the class of extreme events which we are interested in. The land 
use of the surrounding area of a river is usually determined through a long 
period of (non-scientific) observations. The most natural observable is the 
maximum height of the water in a period (e.g., annual maxima), even though 
observables like the drainage area or the discharge volume are also commonly 
used in the scientific literature. After a long period of normality, i.e., when 
the river does not overcome the standard preventive barriers, protection is 
usually neglected, and the barriers inevitably assume a lower value. On the 
opposite, after the occurrence of a flood (extreme event) a lot of attention 
and efforts are directed to avoid similar catastrophes in the future and the 
barriers thus increase. This is the most natural unplanned human reaction 
to extreme events and will constitute the main motivation for the simplified 
model analyzed in this paper. More subtle reactions may affect the observable 
used to characterize the system. For instance, excavations in the river and the 
construction of new buildings and levees in the floodplains, which also depend 
on the occurrence of recent floods, modify the available area of the river. In 
this case, the measure of the water level is directly affected by these human 
activities and not only by the amount of water in the rivers basin. An even 
more drastic human activity can change the dynamics of the system in the 
phase space: the construction of a water reservoir upstream can directly con- 
trol the level of the waters, or, more indirectly, the precipitation in a region is 
influenced by the presence of strong human activity. 

In summary, the human activities act in a kind of feedback loop with the 
occurrence of extreme events and can in principle influence three different 
levels on the measurement chain: 

(11) The preventive barriers, e.g., by increasing the protections around the river. 

(12) The observable used to characterize the system, e.g., by digging the river. 

(13) The dynamics in the phase space in a more fundamental way, e.g., by con- 
structing water reservoirs. 

The reactions to extreme events may be planned or involuntary and, corre- 
spondingly, the two fundamental questions are: 

(Ql) Which is the best method in order to reduce the number of extreme events? 
(Q2) Which is the influence of the feedback reactions on our perception and on 
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the occurrence of extreme events? 

This paper address questions (Ql) and (Q2) through the analysis of a simple 
stochastic model that simulates the feedback reactions (II) and (12). It is orga- 
nized as follows. In Sec. 2 we present the model. The time between successive 
extreme events is studied in Sec. 3 and the efficiency of the model in Sec. 4. 
In Sec. 5 we discuss how our model adapts to the effect of non-stationarities. 
Finally, we summarize our conclusions in Sec. 6. 



2 Moving threshold model 

We consider here a simplified model for the feedback reactions to extreme 
events that takes into account the main features discussed in the previous sec- 
tion. A random sequence of events £ n is taken as a stochastic input to our model 
and represent the complex phenomenon measured in the physically relevant 
observable. We say that an extreme event occurs at time n when £ n overcomes 
the value of the barrier q n , i.e., £ n > q n . In this case we expect the new value 
of the barrier q n+ \ to be increased proportionally to the extreme value £„. On 
the other hand, if no extreme event occurs in time n, the barriers decrease to 
a fraction of its previous value. This decay of the barriers occurs typically due 
to the short memory underlying the human activities (forgetting), but it can 
also appear naturally, e.g., the decay of immunity after vaccination (0), or the 
increasing vulnerability of forest to wind gusts due to the growth of trees. The 
change of the barrier size q n+ i can be thus summarized as 

max{a£ n , f3q n } if £ n > q n , 

Qn+l = < (1) 

[ (3q n if £n < q n , 

where formally < (3 < 1 and a > 0. The max in the first equation can as- 
sume the value j3q n only for a < (3, and is introduced in the model to avoid the 
artificial reduction of the barrier after an extreme event. We study the tem- 
poral sequence of extreme events as a function of the control parameters a, j3, 
with special interest for the cases /3 J$ 1 and a ^ 1, which means that an 
event of the size of the last extreme should not overcome the threshold in the 
(near) future. In principle the dynamics defined by Eq. (1) can be applied to 
any time series that does not contain the influence of human activities. 
In order to avoid further complications of our model we consider initially {£ n } 
to be a Gaussian delta-correlated random variable with (£„) = 0, cr^ = 1, and 

thus p(£) = ~~ and (£n£m) = 8n,m, where (.) denotes temporal average. 

It is interesting to compare this model to other simple stochastic models used 
to simulate, e.g., the occurrence of earthquakes (0), the spikes in neurons (0), 
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and paradigmatic examples of stochastic resonance (J8|). The novel aspect of 
the model studied in this paper is the existence of a threshold that varies 
deterministically in time depending only on the previous extreme events. 

It is sometimes convenient to study the dynamics of Eq. (1) using the variable 



Vn = in~q n , (2) 

where extreme events occur for y n > 0. The mean value and variance of y n 
can be written as 



(y) = (0 - <<?> = -<<?>, 

, r- (3) 

°* =^(y 2 -(y) 2 ) = ^I + ° 2 q -m) 

= V /I +^ >1 ' 

by noting that the term {^q) is zero due to the lack of correlation between £ n 
and q n . 

We would like at this point to associate explicitely our model with the per- 
spective of floods in rivers mentioned before. Considering the original vari- 
ables (£,<?), we regard the human influence restricted to the delimitation of 
the river domain. In this case £ could be the water level and q the size of 
the preventive barrier, measured as the maximum acceptable height of the 
water before causing damage. On the other hand, if we perform the change 
of variable (2), we interpret y as the departure of the water level from this 
threshold (y < below and y > above threshold), while £ as the water in 
the river basin and q as a measure of the modification of the river shape due 
to human activity. We see thus that the dynamics defined by Eq. (1) models 
simultaneously reactions (II) and (12) mentioned in the introduction. 

A general picture of our model is presented in Fig. 1, where numerical results 
of the time series {y n } arid {q n } are shown for three different control param- 
eters a and (3. For typical values (a, (3 ~ 1), the probability density function 
(PDF) p(q) can be approximated by a Gaussian, what leads to a Gaussian 
form of p(y). In this case the knowledge of (y),(Ty (see also Eq. (3)) uniquely 
determines the fraction of extreme events p(y > 0). Increasing a we notice 
an increase of a q and a y . For large a, the distribution p(q) becomes highly 
asymmetric and a long tail for large q's appears. In this case p(y) also loses 
its normal form. 
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Fig. 1. (Color online.) Analysis of the time series of y n = £ n — q ri and q n for three 
different control parameters. For all graphs: black circles represent the PDF of y, 
solid line the Gaussian distribution (with (y) = — (g), cr = cr y ), and histogram the 
PDF of q n . In the insets, time series of q n (dashed line) and y n (solid line with 
circles) are shown. 



3 Interval between extreme events 



One of the most important characteristics of the temporal sequence of extreme 
events is their recurrence time (Jjj), i.e., the time T between two successive ex- 
treme events quantified by the interevent time distribution P(T). Using a 
constant barrier q n = q the extreme events obtained from an uncorrelated 
random time series occur also completely at random and P(T) decays expo- 
nentially. We show in this section that this is not the case when the size of the 
barrier dynamically changes according to Eq.(l). In this case P{T) presents 
typically a maximum, i.e., there is a characteristic interevent time T max > 0. 

To compute P(T), i.e., the probability of having two consecutive extreme 
events separated by time T, we first have to calculate the probability r(t) of 
having one extreme event at time t independent of the other events. In our 
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model the probability of occurrence of one extreme event > qt at time t is 
given by 



r°° r°° 1 t 2 1 \/2 

r(f)=/ P(0de=/ -= e -Tde=-erfc(^- ft ), (4) 

where erfc(x) is the complementary error function. During the interval between 
extreme events the barrier evolves as qt = The value £~ determines go = 

a(,~ and corresponds to the value of £ n that generated the previous extreme 
event, or, in the exceptional case when > q n but (3q n > a£ n [see Eq. (1)], to 
£~ = (3q n /a. Due to the lack of further correlations between extreme events, 
the interevent time distribution P(T) is obtained as the composition of the 
probability of having an extreme event at time T with the probability that no 
event occurred for t = [0, T[, which can be written as 

P(T) = r(T) nV " r(s)) « r(T) exp[]T r(s)), 

S=l 8=1 

where we have used the approximation of small extreme event probability r(T) <C 
1. In the limit of continuous time we obtain (fioh 

P(T) = Cr(T)e~ r{s)ds } (5) 

where C is a normalization constant. Introducing the expression (4) in (5) we 
obtain 

P/T r , C y/2 T T V2a^P T 113 3 1 ^ 2ft iT^ 

P(T-i ) = -erfc(-< )eM-2-^^2F 2 (-,- 2 ,-,----(at )fi )], 

(6) 

where 2P2 () is the hypergeometric function with parameters a\ = a% = \ 
and b\ — b% — | The interevent time distribution for given parame- 

ters a, P, is given by P{T) = J °° P(T\ £~)p{£~)d£~ . However, the distribution 
p(£~) is unknown. To obtain simplified theoretical curves we have inserted in 
Eq. (6) the constant value £ _ = obtained numerically. These distribu- 

tions are plotted in Fig. 2 where we notice the existence of a nontrivial most 
probable interevent time T max , in good agreement with the numerical results. 
This constitutes the main and at first sight most striking result, i.e., extreme 
events occur almost periodically if a > 1. 

The existence of such a most probable interevent interval T max resembles re- 
sults obtained in models presenting stochastic resonance (JSj) or coherence res- 
onance However, this is not the case of our model since it has neither a 
periodic input signal nor a resonance behavior for different noise amplitudes. 
In fact, in our case T max varies with the control parameters a, (3, whereas a 
modified (constant in time) variance of £ n is equivalent to a simple rescaling of 
the length scale and do not affect time scales. We would like to obtain now the 
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Fig. 2. (Color online.) Interevent time distribution P{T) for (3 = 0.95 
and a = 0.8, 1.2, and 4. Lines with symbols are the numerical results and sim- 
ple lines are the theoretical distribution (6) using £~ = (£) = 1.17, 1.66, 1.765, 
respectively. After the maximum at T max , P{T) decays faster than exponentially 
to zero (only times T > have physical meaning). The inset shows the PDF of £~ 
obtained numerically for the same parameters. 

dependence of T max on a, (3. A direct analysis through the theoretical distri- 
bution (6) is difficult due to its complicated expression and due to the lack of 
knowledge of the distribution of £~, which also strongly depends on a, j3. For- 
tunately some intuition can be gained through simple approaches developed 
in what follows. Qualitatively, we notice that the existence of a most probable 
interevent time is a consequence of the reduction of the probability of short in- 
terevent times due to the increment of the barrier size after one extreme event. 
We expect thus that T max increases with a and that for a < 1, P(T) decays 
monotonically with T. These results are verified numerically in Fig. 3a,b. It is 
also interesting to note that the characteristic interevent time T max also shows 
up in the spectrum and autocorrelation function of the series {y n } and {<?„}. 

Consider now that £ n assumes the constant value £* > 0. In this simple case the 
time between events is given by the time the barrier takes to decay to a value 
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Fig. 3. (Color online.) Numerically obtained interevent time distribution P(T) for 
(a) a = 1.2 and p = 0.7,0.85,0.89,0.94,0.98 (from left to right), and (b) for = 0.95 
and a = 0.8, 1, 1.2, 2,4, 10 (from left to right), (c) Auto-correlation function of the 
series y n . (d) Dependence of T max on the control parameters a,/3. 



smaller than £*: q* = £* = a£*f3 T =>- T = — j^rf- Another simple approach is to 
consider that the next extreme event occurs when the barrier is at q + and the 
previous one occurred due to a value £~, where the distributions of (£~,<2 ,+ ) 
are unknown and depend on a, /3. We obtain in this case an interevent time T, 



T 



In a ln(g + ) — ln(£ 



In/? 



ln(/3) 



— rr^r- In both cases we see that 

ln(j9) 

In Fig. 3d we plot numerical 



ln(/3) ' 



On average £~ > q + and thus T max > 

Tmax depends explicitely on the ratio - 

obtained values of T max against — jttj- We notice that all points collapse 
approximately in a same curve that is indeed always above the diagonal and 
that good agreement is achieved by a linear fitting. 
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4 Efficiency of the model 



The most relevant issue in the socio-economic context is how to minimize the 
number of extreme events, or how to optimize the efficiency of a protection 
strategy. The rate of extreme events p(y) > is the simplest and most natural 
measure of the costs due to extreme events and is thus used throughout this 
paper. For many situations it is a realistic measurement of the damage that 
occur due to the overpass of a threshold (e.g, obstruction of a street or of 
the power supply), where there is no difference between "small" and "large" 
extreme events. In other situations, more detailed accounts of the costs would 
make use of cost functions that depend (non-linearly) on y. Our model can be 
considered as a specific method which can be optimized by the choice of the 
control parameters a, (3. 

The rate of extreme events p(y > 0) as a function of the control parameters 
is shown in Fig. 4. Since the barrier q do not necessarily increase after one 
extreme event if a < 1, we notice again a qualitatively different behavior 
for a < 1 and a > 1. For a fixed (3 and varying a (Fig. 4a) we notice that 
the number of extreme events decays drastically around a ~ 1. On the other 
hand, by fixing a and varying f3 (Fig. 4b) the number of extreme events goes 
much faster to zero when (3 — > 1 if a > 1. 

It is quite natural that the number of extreme events is reduced when the con- 
trol parameters increase. However, in real situations the increment of these 
parameters, or equivalently the increment of the barrier, is related to some 
costs that have to be taken into account when studying the efficiency of the 
model. Since the costs are usually increasing with the size of the barrier, we 
measure them by the mean value of the barriers (q n ). In Fig. 5 the rate of 
extreme events is shown against (q) for different values of the control param- 
eters a, (3 and is compared with the result (dotted line) obtained when the 
barrier is maintained unchanged in time. As already suggested in Eq. (3), we 
notice that our moving threshold method is always less efficient than the con- 
stant barrier case. The limit of unchanged barrier is obtained in our model 
for a — > and (3 — > 1, when the most efficient results are obtained. By noting 
that the lines of constant a are approximately parallel to this limit, we realize 
that the relevant limit is f3 — > 1. Indeed, for any a ~ 1 an efficient reduction 
of the number of extreme events is only possible by increasing (3 towards one. 



5 Non-stationarities 

In the previous section we have seen that the model proposed in this paper 
to simulate the feedback reactions to extreme events is always less efficient 
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Fig. 4. Extreme events rate p(y > 0) as a function of the control parameters a and 

than maintaining the size of the barrier constant in time, i.e., a non-reactive 
model. On the other hand, a clear advantage of reactive models is their ability 
to deal with non-stationarities in the time series. This is a specially important 
issue when considering extreme events since in many cases they are indeed 
originated from process presenting slow trends. Once more, these trends may 
be natural or consequence of human activities (e.g., change in land use, global 
warming) . 

In order to explore how the model defined by Eq. (1) adapts to weak non- 
stationarities, we choose in this section the input time series {£„} to be a 
Gaussian delta correlated random variable with mean and variance changing 
linearly in time 



(0 = <£n, 



4-i 



(7) 

a t = 1 + ^T n ' 

where N is the total observation time and (£)^,£r| > are constants. Note 
that increasing the value of (£) is not equivalent to a simple translation since 
the barrier is still limited to positive values q > and lim^oo f3q n = 0. In 
Fig. 6a,b we show that on large time scales the value of the barrier q n also 
increases linearly in time in both cases, i.e., when the mean or the variance 
increases in time. The linear increment of (q) increase also the fluctuations 
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<q> 

Fig. 5. (Color online.) Extreme events rate p(y > 0) as a function of the mean 
barrier (q) (costs) . Dashed lines are obtained by fixing the control parameter j3 and 
varying a, and continuous lines by fixing a and varying [3. The dotted line at left 
correspond to a constant barrier in time qt = (q). 

(o"q) of the time series {q n } and {y n }- 

It is also interesting to compare the results for the interevent time distribution 
of non-stationary time series with those reported in Sec. 3 for the stationary 
case. When (£ n ) increases (decreases) we note that the value T max of the peak 
of the interevent time distribution P(T) slightly decreases (increases). This 
effect is quite natural since the probability of a large value of £ is constantly 
increasing (decreasing) when (£) increases (decreases). On the other hand, due 
to the linearity of Eqs. (1), a change of the variance lead to a rescale of q 
without changing the value of T max . Both effects are verified numerically in 
Fig. (6)c. 

Since the barrier increases proportionally to the size of the extreme event, we 
see that the time our model takes to adjust to non-stationarities is given by 
the interevent time T. When T max is much smaller than the total observation 
time, as considered here, the non-stationary effects can be considered small 
during this time interval. As a consequence of this fast adaptability of our 
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Fig. 6. (Color online.) Time series of the barrier q for a = 1.2,(3 = 0.95 and: (a) 
linear increment in time of (£) with = 5 and (b) linear increment in time 

of o"£ with aj; = 5. The horizontal lines correspond to (q) at n = (bottom) and 
at n = N = 2 15 (top). The insets show magnifications, (c) Numerically obtained 
interevent time distributions for a = 2, (3 = 0.95 and different non-stationary sce- 
narios (see legend). Vertical lines are located at the value of T max . The stationary 
case for a = 2, f3 = 0.95 is indistinguishable from the case where only <7g increases. 



model to the application of relations (7) we have that the dependence of the 
number of extreme events on the parameters a, (3 is qualitatively equivalent 
to the one reported in Sec. 4 for the stationary case. In order to obtain the 
efficiency we have to compare again the extreme events rate p(y > 0) with 
the mean barrier (q) (costs). However, now the value of q is driven by changes 
of (£) (the value (q) reflects only the period of larger (£)) showing that the 
efficiency analysis does not make sense in this case. More interesting is the case 
when the variance changes in time and the mean value (£) is kept constant, 
shown in Fig. 7 for both increasing and decreasing a^. The comparison with 
the results obtained with constant barrier shows that with reasonable choice 
of parameters a, (3 the moving threshold model leads to a much more efficient 
result. 
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Fig. 7. (Color online.) Extreme events rate p(y > 0) as a function of the mean bar- 
rier (q). Dashed lines are obtained by fixing the control parameter /3 and varying a, 
and continuous lines by fixing a and varying 0. The dotted line correspond to a 
constant barrier in time qt = {q). Two non-stationary situation are considered (a) 
<7jf = 5 and (b) = 0. (£) = is constant for both cases. 

6 Conclusion 

In summary, we have introduced a model to simulate human reactions or natu- 
ral feedback response to extreme events. We have obtained that the sequence of 
extreme events occur with a certain periodicity, exclusively due to the human 
activity (natural feedback). Regarding the efficiency of the model, we have 
obtained that the best strategy in order to efficiently reduce the number of 
extreme events is to try to avoid the decrease of the protection barriers in the 
periods between extreme events. On the other hand, if slow non-stationarities 
are present in the phenomena, it is also useful to increase the usual protec- 
tions to the value of the previous extreme event. The same conclusion is also 
expected for positively correlated sequences of events. 

These results are obtained in a very simplified model that tries to isolate the 
influence of the human reactions to extreme events. In more realistic setups the 
properties discussed here may appear together with system-specific character- 
istics. In this sense, we can relate the characteristic interevent time observed 
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in our model with the observed interepidemic interval between, e.g., small- 
pox epidemics (flih. Direct association of realistic preventive schemes with the 
control parameters of our model lead to the estimation of the characteristic 
interevent time. For instance, in the example of floods in river we may esti- 
mate that barriers are reduced by 2% every year and that after a flood they 
are increased by 20% more than the highest water level. With these parame- 
ters and ignoring deviations of Gaussianity, we obtain through our model the 
reasonable estimation of a 53 years period between floods. 
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